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MOND and the Lensing Fundamental Plane: No need for 
dark matter on galaxy scales 
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ABSTRACT 

Bolton ct al. (2007) have derived a mass-based fundamental plane using photometric 
and spectroscopic observations of 36 strong gravitational lenses. The lensing allows 
a direct determination of the mass-surface density and so avoids the usual depen- 
dence on mass-to-light ratio. We consider this same sample in the context of modified 
Newtonian dynamics (MOND) and demonstrate that the observed mass-based fun- 
damental plane coincides with the MOND fundamental plane determined previously 
for a set of polytropic spheres chosen to match the observed range of effective radii 
and velocity dispersions in elliptical galaxies. Moreover, the observed projected mass 
within one-half an effective radius is consistent with the mass in visible stars plus 
a small additional component of "phantom dark matter" resulting from the MOND 
contribution to photon deflection. 



1 INTRODUCTION 



Modified Newtonian dynamics (MOND) posits a fundamen- 
tal acceleration scale (ao) below which the law of gravity or 
inertia deviates from the Newtonian form (Milgrom 1983). 
Therefore, a generic prediction is that in systems where the 
internal acceleration, expressed as a surface density, exceeds 
a critical limit (ao/G), there should be little discrepancy 
between the observable mass and the Newtonian dynamical 
mass. In traditional language, there should be no evidence 
for dark matter within high surface brightness systems. 

Examples of high surface brightness systems are globu- 
lar clusters and luminous elliptical galaxies. It is well-known 
that the standard Newtonian analysis of the internal kine- 
matics of globular clusters yields mass-to-light ratios that 
are completely consistent with normal stellar populations 
(e.g., Pryor et al. 1989). However, dark matter proponents 
were surprised by the more recent analysis of planetary neb- 
ulae kinematics in several luminous ellipticals which implied 
a "dearth of dark matter" within the bright inner regions of 
these galactic systems (Romanowsky et al. 2003). This re- 
sult is completely consistent with the expectations of MOND 
(Milgrom & Sanders 2003). 

The galaxy scale gravitational lenses that produce mul- 
tiple images of background sources- strong lenses- are pri- 
marily in early type galaxies: ellipticals or bulge dominated 
spirals (Kochanek et al. 2004); that is to say, high surface 
brightness systems. Because the surface density of lenses 
necessary to produce multiple images generally exceeds the 
MOND critical surface density, a MOND corollary predic- 
tion is that the implied mass surface density in strong grav- 
itational lenses should be generally consistent with that of 
observable starlight, i.e., little or no dark matter in strong, 
galaxy scale, gravitational lenses, at least within the Ein- 
stein ring radius: "what you see is all there is." 



Recently, Bolton et al. (2007) have considered the "fun- 
damental plane" of early type galaxies (the observed rela- 
tion between effective radius, surface brightness and velocity 
dispersion) as defined by a sample of 36 strong lens galaxies 
from the Sloan Lens ACS. The lensing analysis was com- 
bined with spectroscopic and photometric observations of 
the individual lens galaxies in order to generate a "more 
fundamental plane" based upon mass surface density rather 
than surface brightness. They found that this mass-based 
fundamental plane relation exhibits less scatter and is closer 
to the expectations of the Newtonian virial relation than 
is the usual luminosity-based fundamental plane. This pre- 
sumably is due to the absence of both an intrinsic scatter 
in M/L as well as its systematic variation with luminosity. 
Below, we emphasise that the implied lensing M/L values 
within the Einstein ring radius do not require the presence 
of a substantial component of dark matter, again consistent 
with the expectations of MOND. 

With MOND, the mass-velocity dispersion (M — a) re- 
lation (or its observational equivalent, the luminosity-based 
Faber- Jackson law), is, in some sense, more fundamental 
than the fundamental plane. In fact, near isothermal, ho- 
mologous objects exhibit a precise mass-velocity dispersion 
relation (Milgrom 1984). But, in order to match the observed 
range of elliptical galaxy properties- primarily the distribu- 
tion of the velocity dispersion vs. effective radius- models 
for such systems must deviate from strict homology. The 
necessary deviation yields a large scatter in the predicted 
M — a law; none-the-less, the models define a narrow mass- 
based fundamental plane when an additional parameter is 
included, such as effective radius or surface density (Sanders 
2000). 

Here, we compare this previously derived MOND funda- 
mental plane to the mass-based fundamental plane of Bolton 
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Figure 1. The points illustrate the mass- based fundamental 
plane of Bolton et al. (2007) shown as a relation between the 
observed projected mass within r e f t /2 and the dimensional mass 
(eq. 2). Also shown is the same for the MOND fundamental plane 
as determined from 360 anisotropic high n polytropes (Sanders 
2000). Mass is given in units of 10 11 M Q . 

fit by 

log(M iens ) = 5\og(M dim ) + C. (3) 

with 5 — 0.98 and C « 0.6. Therefore, these observations 
are consistent with the Newtonian virial theorem and no 
variation of the structure constant; elliptical galaxies would 
appear to be quite homologous over a wide range of mass. 
Moreover, the lensing M/L values (r-band) for these galaxies 
range from about two to eight in solar units and are some- 
what higher than predicted by population synthesis models 
(Fig. 3, upper panel). This point is addressed further below, 
but, overall, for early-type galaxies M/L values in this range 
would not constitute compelling evidence for dark matter 
within r e ///2. 



3 THE MOND FUNDAMENTAL PLANE 

By solving the structure equation, Milgrom (1984) has 
demonstrated that, with MOND, isothermal spheres with 
a fixed degree of anisotropy have finite mass, an upper limit 
to the surface density (« ao/G) and exhibit a M — a 4 re- 
lationship (the basis of the Faber- Jackson law); therefore, 
such objects might constitute sensible models for ellipti- 
cal galaxies. However, the observed properties of ellipticals- 
in particular, their distribution and scatter on the cr-r e ff 
plane (J0rgcnscn 1999. J0rgcnsen, Franx & Kargard 1995a| 



et al. and find that they are entirely consistent. Using the 
theoretical fundamental plane relation to calculate the mass 
from the observed effective radius and velocity dispersion, it 
is found that the projected MOND FP mass, presumably en- 
tirely baryonic, is proportional to the observed lensing mass 
with a small offset due to the effect of modified gravity on 
the deflection of photons along the line of sight. Moreover, 
the implied MOND M/L-colour relationship for these sys- 
tems is entirely consistent with population synthesis models. 
In other words, we find that MOND requires no dark mat- 
ter on galaxy scale in lensing galaxies, in contrast to recent 
claims (Ferreras et al. 2008). 



2 THE MORE FUNDAMENTAL PLANE 

The fundamental plane of elliptical galaxies (Dressier et al. 
1987, Djorgovski & Davis 1987) is a scaling relationship in- 
volving the observed effective radius, r e //, the central line- 
of-sight velocity dispersion, ao, and, in the original form, the 
mean surface brightness within an effective radius, I. This 
is usually expressed as 

log(r e ) = a log(cr) + b log (J) + d (1) 

(it may alternatively be written as a relationship between 
total luminosity, effective radius, and velocity dispersion). 
The analysis of Sloan data on elliptical galaxies implies that 
a w 1.5 and b w —0.8 (Benardi et al. 2003), while the expec- 
tations from the Newtonian virial theorem, assuming ho- 
mology and constant M/L, would be a = 2 and b — — 1. 
The observed deviation from these expectations is generally 
thought to be due to systematic deviations from homology 
or from constant M/L. 

The significant achievement of Bolton et al. was to use 
gravitational lensing in order to eliminate the dependence on 
M/L. The strong gravitational lensing provides the mass sur- 
face density within the Einstein ring radius (r e ; n ); therefore 
surface density, E, may be substituted for surface brightness 
in the above relationship. The Einstein ring radius in their 
sample is generally about half the effective radius, so the 
surface density is taken by Bolton et al. to be that within 
r e ///2 (the correction from the r e i n to r e ///2 is done either 
by assuming a mass distribution like that of an isothermal 
sphere or that light traces mass). This was then combined 
with spectroscopic and photometric observations of the lens 
galaxies to generate a "mass" fundamental plane. The re- 
sulting relation exhibits considerably less scatter and the ex- 
ponents are close to the virial expectations: a — 1.86±0.17, 
b = -0.93 ± 0.09, d = 5.4 ± 0.9 with the mass traces light 
assumption (here the effective radius is in kpc, the velocity 
dispersion in kms -1 and the surface density in M Q kpc -2 ). 

The Newtonian virial theorem implies that the mass 
within radius r is given by 

M = ca 2 r/G (2) 

where c is a constant determined by the structure of the 
object. Therefore Bolton et al. also express their result by 
plotting the mass of the lens (within r e ///2) against a "di- 
mensional" mass variable given by a 2 r e ///(2G) where a is 
the mean line-of-sight velocity dispersion within r e ///2. The 
result of this is shown by the points in Fig. 1 which are well 



|j0rgensen, Franx fc Kaerg ard 1995bJ- cannot be matched 
by isothermal spheres or any strictly homologous class of 
objects. The MOND isothermal sphere is too inflated for a 
given velocity dispersion- the mean surface density within 
an effective radius is too low. In order to match these ob- 
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served properties of elliptical galaxies with MOND, it is nec- 
essary to exploit other degrees of freedom: models must de- 
viate from being strictly isothermal with a fixed degree of 
orbital anisotropy. 

In order to explore these possible degrees of freedom, 
Sanders (2000) has considered polytropic spheres with an 
anisotropy parameter, /3 = 1 — a-t 2 /a r 2 , which varies sys- 
tematically with radius as in the Osipkov-Merritt models, 
(3 — r 2 /(r 2 + r a 2 ), (r a is the anisotropy radius, Binney & 
Tremaine 1987) . For a polytropic sphere of index n the radial 
velocity dispersion-density relationship is given by 

*r 2 = Ap 1/n (4) 

where A is a constant (n — *■ oo corresponds to the isother- 
mal sphere). It was found that a range 12 < n < 16 was 
sufficient to match the the mean value and wide disper- 
sion of effective radius for a given velocity dispersion, pro- 
vided that the anisotropy radius also varies over the range 
2 < r a /r e ff < 25. The luminosity or mass density distribu- 
tion within these objects is similar to that of a Jaffe model 
(Jaffe 1982): within roughly r e // the density falls as r~ 2 
and beyond steepens to r~ 4 . Spherically symmetric N-body 
calculations with MONDian modified gravity (Bekenstein 
& Milgrom 1984) demonstrate that objects resembling such 
large n polytropic spheres may actually condense and recol- 
lapse out of the Hubble flow (Sanders 2008). 

Each model, characterised by a particular value of n and 
r a /r e ff, exhibits its own exact M oc a 4 relationship but for 
all models combined there is a large dispersion in this rela- 
tion. None-the-less, in spite of the dispersion in homology, 
the models lie on a narrow fundamental plane 

log(r e ) = 2.0 log(a e2 ) - 1.06 log(E e2 ) + 6.16 (5) 

where <r e 2 and £ e 2 are the velocity dispersion and the av- 
erage surface density respectively within r e ///2 (the coeffi- 
cients here differ slightly from those given by Sanders be- 
cause of scaling to r e ///2). This is seen to be quite close to 
the Newtonian virial expectation and within the errors of the 
more fundamental plane defined by the gravitational lenses. 
In fact, both the MOND and lensing fundamental planes are 
well-approximated by the virial relation for Newtonian Jaffe 
models. 

As above, the MOND fundamental plane relation may 
be expressed in terms of a MOND FP mass (projected 
within r e ff/2) as a function of the dimensional mass unit 
a 2 r e ff /2G. This is shown by the thick line in Fig. 1. This is 
the original ensemble of 360 anisotropic, polytropic models 
covering the range described above. We see that it is almost 
coincident with the observations of Bolton et al. 

We may also apply the MOND fundamental plane rela- 
tion in order to calculate the mass of a lens galaxy, projected 
within r e ///2, from the observed velocity dispersion and ef- 
fective radius. This calculated MOND FP mass is plotted 
against the observed lensing mass in Fig. 2. It is evident 
that the MOND FP mass is closely proportional to the lens 
mass, but about 30% less on average. This is understand- 
able because the MOND FP mass is the total projected 
baryonic mass; the lensing mass, however, includes a phan- 
tom dark matter contribution because MOND, or its rela- 
tivistic extension, TeVeS (Bekenstein 2004), provides extra 
deflection along the line-of-sight (Mortlock & Turner 2001 
IZhao et al. 2006(1 . 
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Figure 2. The logarithm of the observed lensing mass at r e j j /2 
(Bolton et al. 2007) plotted against the logarithm of the projected 
MOND FP mass within r eff /2 calculated via the MOND funda- 
mental plane, eq. 5 (Sanders 2000). The solid line is the line of 
equality. The 30% offset is obvious. 
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Figure 3. The upper panel shows the observed lensing M/L 
within r e ff/2 in the r band plotted against the r — i colour. 
The solid line is the prediction of the population synthesis mod- 
els of Bell et al. 2004 with the Kennecut or Kroupa initial mass 
function. The dashed lines show the 0.2 dex range which is a re- 
alistic estimate in the uncertainty of the estimated stellar M/L 
values given the uncertainty of the initial mass function and the 
effects of metallicity. The lower panel is the same for the MOND 
M/L. 
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It is of interest to compare the MOND M/L values with 
population synthesis models. In Fig. 3 (lower panel) the in- 
dicated M/L values within r e ///2 are plotted against r — i 
colour and compared to the theoretical models of Bell et 
al. (2003). The luminosities have been determined from the 
observed Sloan magnitudes applying both K- and evolution- 
ary corrections (Poggianti 1997) and reduced by a factor of 
about three to correspond to that within r e ///2. We see 
that MOND M/L values closely track the theoretical val- 
ues. The upper panel is the same except that the M/L is 
that determined from the projected lensing mass. Again the 
same trend is present although the offset discussed above is 
evident. 



4 THE CONTRIBUTION OF MODIFIED 
GRAVITY TO THE DEFLECTION OF 
PHOTONS 

The current multi-field relativistic extensions of MOND, 
such as TeVeS, are characterised by a physical metric which 
is distinct from the Einstein, or gravitational, metric. The 
transformation between the two metrics is "disformal" and 
chosen such that the relationship between the total weak 
field gravitational force and the deflection of photons is iden- 
tical to that of General Relativity; specifically, the deflection 
angle is given by 



g±dz 



(6) 



where the integral is along the line of sight and g± is the 
perpendicular component of gravitational force. The grav- 
itational force includes not only the traditional Newtonian 
force but also an additional 'MOND' force (mediated by a 
scalar field) which becomes dominant at low accelerations. 

This MOND force may be viewed as arising from a 
phantom dark halo. The phantom halo begins to make a 
significant contribution to the gravitational force and hence 
to the deflection of photons near the MOND transition ra- 
dius 



yj GM / ao 



(7) 



where M is the true (baryonic) mass of the object. In fact, 
the phantom halo may exhibit an apparent shell of matter 
between 1/2 r t up to rt depending upon the form of the 
function which interpolates between the MOND regime and 
the Newtonian regime (Milgrom & Sanders 2008). The rel- 
evance for the present discussion is that we would expect to 
find some evidence of "dark matter" in strong gravitational 
lensing if the Einstein ring radius is a significant fraction of 
the MOND transition radius; the larger that fraction, the 
larger the apparent contribution of dark matter to the pro- 
jected mass detected by strong gravitational lensing. 

The Einstein ring radius, for a point mass, is given by 



4GM DiDu 
c 2 D s 



(8) 



where Di, D s , and Di s are respectively the angular size dis- 
tances from the observer to the lens, to the source and from 
the lens to the source. Combining eqs. 7 and 8 and making 
use of the fact that ao = fcHo where / w 1/6 we find that 
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(9) 
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Figure 4. The Einstein radius in terms of the MOND transition 
radius as a function of the lens redshift for three different values 
of the source redshift: z s = 0.5 (dashed curve), 3.5 (dotted), 100 
(solid). This is for the standard "concordance" cosmology. 



where F is a function of the lens and source redshift and 
depends upon the cosmological model. 

This function F(zi,z s ) is shown in Fig. 4 as a func- 
tion of zi for three different values of the source redshift 
(z s = 0.5, 3.5, 100). This is all done for the standard "con- 
cordance" cosmology which the relativistic MOND cosmol- 
ogy should mimic to a close approximation. For the lenses 
in the sample of Bolton et al., z s « 0.5 and zi « 0.2 on aver- 
age; therefore, F does not exceed about 0.25, and we might 
expect the contribution of an apparent halo to lensing to be 
relatively small. On the other hand for strong lenses con- 
sidered by Treu & Koopmans (2004), Zi « 1 and z s w 3.5 
implying F « 0.4. The expectation is that there should be 
evidence for more "dark matter" in these high redshift sys- 
tems, as, in fact, is reported by Treu & Koopmans. 

We have estimated the contribution of modified gravity 
to lensing, in the form of projected phantom dark matter, 
for all of the objects in the sample of Bolton et al. The 
objects have been modelled as isotropic Jaffe spheres with 
the observed effective radius. Given this density distribu- 
tion, the Jeans equation is solved for the run of radial ve- 
locity dispersion where the total force, g is related to the 
Newtonian force, gjv, by the MOND formula gfi(g) = pjv; 
n{x) is the MOND interpolating function (/J.(x) = x where 
x « 1 and /J.(x) — 1 where x >> 1), taken here to be 
of the form suggested by Zhao & Famaey (2006). In each 
case we adjust the mass of the sphere in order to reproduce 
the observed luminosity-weighted average line-of-sight ve- 
locity within r e ///2; this mass in most cases agrees closely 
with that implied by the MOND fundamental plane (eq. 5). 
The difference between total (MOND) acceleration and the 
Newtonian acceleration (\g — (?jv|) then allows us estimate 
the density distribution of the phantom dark halo. The pro- 
jected MOND FP mass (X-axis in Fig. 2), presumably the 
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Figure 5. The observed lensing mass within r e ff/2 from Bolton 
ct al. plotted against the projected MOND FP mass including 
the contribution of phantom dark matter (Mpjjm- This corrects 
Fig. 2 for the effect of modified gravity on gravitational lensing. 

visible mass of the galaxy, is then "corrected" by adding in 
the projected phantom dark mass. The result is shown in 
Fig 5 which again shows the lensing mass vs. the MOND 
FP mass including the phantom dark matter. We see that a 
lensing mass which is 30% higher than the MOND mass is 
explained by the contribution of modified gravity to photon 
deflection. 

We should note, however, that the Zhao-Famaey inter- 
polating function favours the appearance of phantom dark 
mass within the optical image of the galaxy. This is because 
the transition from Newton to MOND is rather more grad- 
ual than for the often assumed ('standard') form of fx applied 
to calculation of galaxy rotation curves (fi(x) = x/\/l + x 2 ). 
Applying the standard form would result in a 10% reduction 
in the projected phantom dark mass, so the appearance of 
Fig. 5, and the conclusions we draw from it, would not be 
altered. 



5 CONCLUSIONS 

With MOND, the baryonic mass-rotation velocity relation 
for spiral galaxies, which forms the basis of the Tully-Fisher 
law, is exact in so far as it relates to the asymptotic rota- 
tion velocity measured far from the luminous galaxy. On 
the other hand, the mass-velocity dispersion relation for 
pressure supported systems, the basis of the Faber-Jackson 
law, is only exact for homologous models; the scaling of 
the relation depends upon the detailed characteristics of 
the object. Actual elliptical galaxies exhibit a range of 
properties- various shapes, varying degrees of deviation from 
an isothermal state and, no doubt, isotropy of the velocity 
dispersion- and cannot be represented by a single homolo- 
gous sequence of models. Therefore, spheroidal galaxies will 



inevitably present a Faber-Jackson law with considerable 
scatter. None-the-less, MOND provides an explanation for 
the remarkable fact that self-gravitating, pressure-supported 
quasi-isothermal objects with a velocity dispersion of a few 
hundred kms -1 will have a mass in the range of galaxies- or 
objects with a velocity dispersion < 10 kms -1 will have the 
mass of globular clusters- or objects with 1000 kms -1 will 
have the mass of a cluster of galaxies. 

In spite of the scatter in the mass-velocity dispersion 
relation, when an additional parameter is added, such as 
effective radius or surface brightness, MOND models for el- 
liptical galaxies define a narrow fundamental plane which 
is close to that implied by the Newtonian virial relation for 
homologous objects (isotropic Jaffe models). This was not 
part of the original set of MOND predictions but became ap- 
parent when it was realized that normal elliptical galaxies 
are essentially Newtonian systems within the effective radius 
and exhibit a wide dispersion in the effective radius- velocity 
dispersion relation. The properties of this fundamental plane 
were outlined by a set of 360 large n polytropic spheres with 
radially dependent anisotropy chosen to match the observed 
joint distribution of ellipticals by effective radius and veloc- 
ity dispersion (Sanders 2000). Applying this fundamental 
plane relation to determine the mass of those ellipticals in 
the sample of J0rgenson et al. (1995) yielded reasonable val- 
ues for the mass-to-light ratios. 

Now, thanks to the work of Bolton et al. (2007) we can 
compare this mass-based MOND fundamental plane directly 
to the observed mass-based fundamental plane as denned by 
this set of 36 strong gravitational lenses. Figs. 1 and 2 illus- 
trate that the two coincide apart from a systematic offset of 
about 30%. Indeed, the implied MOND mass-to-light ratios 
are completely consistent with population synthesis models 
(Fig. 3) , and the small discrepancy between the lensing mass 
and the MOND FP mass can be understood in terms of the 
contribution of modified gravity to the deflection of pho- 
tons (Fig. 5). It is important to recall that the properties of 
the MOND fundamental plane (Sanders 2000) were defined 
well before those of observed mass-based fundamental plane 
(Bolton et al. 2007), so this does, properly speaking, consti- 
tute a prediction that has been subsequently confirmed. 

Most significantly, there is no evidence from strong 
gravitational lensing for a significant mass discrepancy 
within these high surface density systems- as MOND would 
robustly predict. This is in contrast to a recent claim by Fer- 
reras et al. (2008) based upon lensing by six early type galax- 
ies. They note that the lensing mass, as determined either 
by General Relativity or MOND (as extended by TeVeS), 
is significantly greater than the stellar mass estimated via 
population synthesis models. However, this conclusion ap- 
pears to give much weight to the precision of such models; 
the mass difference is generally smaller than the differences 
due to the assumption of different initial mass functions 
(Ri 0.2 — 0.3 dex). Moreover, in the near infrared, the scatter 
induced by metallicity effects can be comparable (Bell et al. 
2003). Overall it is difficult to argue that implied mass-to- 
light ratios ranging from two to eight constitute compelling 
evidence for dark matter in early-type galaxies. 

The conclusions here are generally consistent with those 
of Zhao et al. (2006). Looking at a sample of 18 strong 
double-image lenses from the CASTLES sample in the con- 
text of TeVeS, they concluded that the implied lensing 
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masses are consistent with the observed baryonic matter 
(with additional deflection provided by the modified grav- 
ity). There are two extreme outliers, but this is not unsur- 
prising in view of the possible contribution of surrounding 
groups or clusters to deflection of photons (in one case, the 
implied lensing mass is significantly less than the estimated 
stellar mass which would also be a problem for GR). Treu 
& Koopmans (2004) have argued for the presence of a sig- 
nificant fraction of dark matter (/dm = 0.3 — 0.7) within 
the Einstein ring radius of five well-observed high redshift 
lens systems (for the systems considered in this paper inter- 
preted in terms of dark matter /dm — — 0.4). We have 
argued here that such a result is expected because of the 
larger contribution of modified gravity to the deflection of 
photons in higher redshift systems where the Einstein ring 
radius is a larger fraction of the MOND transition radius 
(Fig. 4). 

The results here for this homogeneous relatively nearby 
sample of lenses are encouraging for MOND- primarily 
with respect to the basic prediction of small discrepancy in 
high surface density systems and secondarily in confirming 
the derived MOND fundamental plane. Such observations 
clearly constitute a powerful test of MOND as an accelera- 
tion based modification. Examples of isolated HSB galaxies 
showing a large discrepancy within the bright regions, would 
be problematic for the theory. 

We thank Leon Koopmans and Adam Bolton for shar- 
ing their data on the SLACS gravitational lenses. We are 
also very grateful to Leon and to Moti Milgrom for helpful 
comments on the initial manuscript. 
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